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Detailed quantitative predictions are obtained for phonon and multiphonon excitations in well- 
deformed rotor nuclei within the geometric framework, by exact numerical diagonalization of the 
Bohr Hamiltonian in an SO (5) basis. Dynamical 7 deformation is found to significantly influence 
the predictions through its coupling to the rotational motion. Basic signatures for the onset of rigid 
triaxial deformation are also obtained. 
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I. INTRODUCTION 

The Bohr Hamiltonian [1, 2], together with its gener- 
alizations [3, 4], has long served as the conceptual bench- 
mark for interpreting quadrupole collective dynamics in 
nuclei. The conventional approach to numerical diago- 
nalization of the Bohr Hamiltonian, in a five-dimensional 
oscillator basis [4-6], is slowly convergent and requires a 
large number of basis states to describe a general de- 
formed rotor-vibrator nucleus. Therefore, it has com- 
monly been necessary to apply varying degrees of approx- 
imation in addressing the dynamics of transitional and 
deformed nuclei, as in the rotation-vibration model [7] 
and rigid triaxial rotor [8] treatments of the Bohr Hamil- 
tonian, or in more recent studies of critical phenom- 
ena [9-12]. 

However, diagonalization of the Bohr Hamiltonian is 
now possible [13] for potentials of essentially arbitrary 
stiffness. In particular, the algebraic collective model 
(ACM) [14-18] provides an efficient and straightforward 
computational framework based on SU(1, 1) x S0(5) al- 
gebraic methods. The Bohr Hamiltonian is diagonalized 
in a basis of SU(1, 1) x S0(5) product wave functions 
on the Bohr deformation variables /3 and 7 and Euler 
angles fi. These are of the form R^{a; f5)'^^aLM{l^^), 
where i?^ is an SU(1,1) modified oscillator wave func- 
tion [19] and ^vaLM is an S0(5) D S0(3) spherical har- 
monic [20, 21]. The S0(5) D S0(3) formulation may be 
used either simply to extend the conventional oscillator 
basis to higher phonon numbers sufficient to provide full 
convergence [22-24] or, further, to obtain much faster 
convergence as a function of basis size through the use 
of SU(1, 1) /3 wave functions chosen optimally for the nu- 
clear deformation [15]. 

The Bohr Hamiltonian can consequently be applied, 
without approximation, to the full range of nuclear 
quadrupole rotational- vibrational structure, from spher- 
ical oscillator to axial rotor to triaxial rotor. Full con- 
vergence can be obtained for energies and electromag- 
netic transition strengths involving high-lying states, for 
instance, intcrband transitions among /3, 7, and multi- 
phonon bands in well-deformed rotor nuclei. The Bohr 
Hamiltonian inherently induces coupling of the /3, 7, and 



rotational degrees of freedom, thereby yielding a rich set 
of phenomena. 

To approach an understanding of the full problem, we 
shall consider, in this article, the simpler but already ex- 
tensive implications of coupling of the 7 and rotational 
degrees of freedom. The relevant Hamiltonian is then 
the "angular" part of the Bohr Hamiltonian, and the 
ACM calculation reduces to diagonalization in a basis 
of S0(5) D S0(3) spherical harmonics (Sec. II). The 
regime we address consists of rotational structure with 
axially symmetric (axial) or weakly triaxial deformation. 
However, even for a nominally axial rotor, the Bohr de- 
scription is found to mandate significant dynamical fluc- 
tuations in 7, far from 7 = 0°. The evolution of spec- 
troscopic quantities (energies and transition matrix el- 
ements) with respect to the 7 confinement provided by 
the potential is systematically investigated (Sec. HI), and 
the spectroscopic implications of the onset of rigid tri- 
axial structure are explored (Sec. IV). Probability dis- 
tributions with respect to 7 and with respect to the K 
quantum number are then used to examine the degree of 
adiabaticity, or separation of rotational and vibrational 
degrees of freedom in the wave functions (Sec. V). Pre- 
liminary results were presented in Refs. [25, 26]. 



II. HAMILTONIAN AND SOLUTION METHOD 



Hamiltonian 



The Bohr Hamiltonian [2] is given, in terms of the 
quadrupole deformation variables j3 and 7 and Euler an- 
gles fl, by 
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The operator appearing in brackets in the kinetic energy 
is the Laplacian in five dimensions. Its angular part A^ 
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is the Casimir operator for the five-dimensional rotation 
group SO (5), which contains the rotations in physical 
space, acting on the Euler angle coordinates, as an S0(3) 
subgroup. The Bohr coordinates are five-dimensional 
spherical polar coordinates, in terms of which the five 
components qm {M = —2, . . ., 2) of the quadrupole de- 
formation tensor are expressed as 



QM = 13 



cos 7, 



-^sin7[^^'i(0) + ^i^)^(n)] 



(3) 



The potential energy V{(3, 7) must be periodic in 7, with 
period 120°, and it must be symmetric about 7 = 0° and 
7 = 60°. The Bohr coordinate system and Hamiltonian 
are reviewed in detail in, e.g., Ref. [27]. 

The restriction to angular coordinates (7, f2) then 
yields a Hamiltonian 



H = A^ + V{j). 



(4) 



Such an angular Hamiltonian arises as a schematic limit 
of the full Bohr Hamiltonian when the coordinate ^5 in (1) 
is taken to be rigidly fixed, as might be considered for 
a well-deformed nucleus. However, a reduction to the 
angular form (4) is more broadly applicable to transi- 
tional nuclei as well [11, 12], since it occurs by sepa- 
ration of variables when the potential is of the form 
1/(^,7) = u{l3) +v{'y)/l5'^ [28]. The explicit relations 
for reduction to an angular Hamiltonian are reviewed in 
Appendix A. The symmetry conditions on 1^(7) are sat- 
isfied by the function cos 87 and powers cos" 87 thereof. 



Let us therefore consider, in particular, 

F = A^ -h x[(l - cos 37) -I- ^cos^ 87] . 



(5) 



The possible shapes of the potential appearing in this 
Hamiltonian are shown in Fig. 1. For ^ = 0, 1^(7) oc 
(1 — cos 87), as considered in Ref. [14], providing a min- 
imum at 7 = 0° (axial deformation). With increasing 
X, a "deeper" potential provides greater confinement or 
stabilization around 7 = 0°, approximately harmonic 
(oc 7^) for small 7. Including a cos^ 87 term [Fig. 1 
(dotted curve)] by taking ^ nonzero introduces a richer 
extremum structure and a means for studying the axial- 
triaxial shape transition [10]. For ^ = 1/2, the potential 
is more softly confining in 7, with a quartic minimum 
(locally cx 7^). This case is termed "critical" in Ref. [10]. 
For ^ > 1/2, the potential has a minimum at a nonzero 
value of 7, given by cos 870 = 1/(2^). For large pos- 
itive ^, the cos^ 87 term dominates, and the minimum 
approaches 7 = 80°. Although not considered here, with 
a negative cos^ 87 contribution the Hamiltonian (5) may 
also be used to investigate prolate-oblate shape coexis- 
tence [29]. 
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FIG. 1: The shape of the potential ^(7) used in (5), plotted 
for various values of £, (taking x = !)■ Note that a constant 
offset 5 has been subtracted from each curve, so that V{Q) = 
in each case. The dotted curve indicates the shape of the 



contribution from cos 87. 



B. Solution method 

Any function of the coordinates (7,^) with the req- 
uisite symmetry properties for a wave function can be 
expressed in terms of symmetric linear combinations of 
Wigner S) functions as {e.g., Ref. [27]) 



V(7,^i)= E^^W4'^m(")' 

even 



(6) 



where [21] 



{1 + 5kY/^ 



(7) 

The wave function is thus fully specified by the F/f (7). 

A complete set for expanding wave functions i/;(7, f2) 
is provided by the SO (5) D SO (8) spherical harmonics 
"^vaLMil,^) [20, 21]. The S0(5) D S0(8) spherical har- 
monics are defined as the eigenfunctions of the SO (5) 
Casimir operator A^, with 



(7, Q) = v{v + 8)*^„lm(7, (8) 



chosen furthermore to posess definite angular momentum 
with respect to the S0(3) subgroup of physical rotations. 
The ^vaLM are labeled by the S0(5) seniority quantum 
number u (t> = 0, 1, . . .), the S0(3) angular momentum 
quantum number L, and its z-projcction quantum num- 
ber M. (A multiplicity index a is also required to com- 
plete the labeling for t; > 6 but will be omitted from the 
notation below when not needed.) The '^vaLM are explic- 
itly realized by constructing the functions Fk (7) needed 
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to express each spherical harmonic in the form (6), as 
may be accomphshed by the algorithm of Refs. [20, 21]. 

Diagonalization of the Hamiltonian (5) is carried out 
in a finite basis of these SO (5) D SO (3) spherical har- 
monics, truncated to some maximum seniority Vmax- In 
general, higher-seniority spherical harmonics are needed 
for the construction of more highly 7-localized wave func- 
tions. Thus, diagonalization for Hamiltonians with stiffer 
7 confinement requires a basis with higher Wmax- A ba- 
sis with z;max = 50 amply suffices for convergence of all 
calculations in the present work. 

It is first necessary to compute the Hamiltonian ma- 
trix elements with respect to the S0(5) D S0(3) 
basis. For the kinetic energy, the matrix elements 
(^'i)'a'LM|A^|\E'„c(LM) are trivially evaluated by the eigen- 
value equation (8). For the potential energy, the ma- 
trix elements of cos 37 may be evaluated in terms of 
integrals of products of Fk{j) functions [14]. Since 
^'30o(7,f^) = (87r2)-i/2(3/y^)cos37, it may be noted 
that the matrix elements of interest are triple overlaps 
{^v'a'LM\^30o\^vaLM) of Spherical harmonics, which are 
equivalent to SO (5) D SO (3) generalized Clebsch-Gordan 
coefficients [20, 21]. These are calculated and tabulated 
electronically (for v < 50) in Ref. [21]. The matrix ele- 
ments of cos" 37 follow immediately from those of cos 37, 
by insertion of resolutions of the identity, i.e., by matrix 
multiplication. 

Then, diagonalization of the Hamiltonian matrix yields 
the amplitudes Oz,ij in the decomposition 

i^LiMil,^) = ^auj'bLjM{l,^)- (9) 

i 

Here we have denoted the zth eigenfunction of the Hamil- 
tonian, for angular momentum L, by 'tpuMi'i and 
likewise relabeled the jth S0(5) D S0(3) spherical har- 
monic of angular momentum L as ^ljm, i-e-, replacing 
V and a by a simple running index [21]. 

The leading-order electric quadrupole operator in 
the Bohr framework is A4{E2) oc q. Under the 
present restriction to angular coordinates, A4{E2) oc Q, 
where Q is the unit quadrupole tensor [20], defined by 
9m = I3Qm [see (3)]. It is straightforward to calcu- 
late transition matrix elements between the Hamilto- 
nian eigenstates (9), once the matrix elements are ob- 
tained between the basis states. Since 'i^i2M{'y,^) = 
(87r^)^-^/^\/15/2QM, the reduced matrix elements are 
proportional to {"^v'a'L'\\'^i2\\'^vaL), which are again 
given by S0(5) D S0(3) generaUzed Clebsch-Gordan co- 
efficients, available from Ref. [21]. 

III. PHONON AND MULTIPHONON 
EXCITATIONS 

A. Spectra 

The nature of the spectra obtained from the Hamil- 
tonian (5) depends both on the depth of the potential 



(determined by x) and the shape of the potential (deter- 
mined by ^ as in Fig. 1). The depth of the potential ef- 
fectively controls the degree of 7 confinement. It is worth 
first carefully considering the implications of 7 confine- 
ment, or conversely 7 softness, within this Bohr Hamil- 
tonian framework. In this section, we shall therefore in- 
vestigate the structural dependence on x (for ^ = 0), 
before proceeding to the dependence of structure on the 
shape of the potential, and in particular the onset of rigid 
triaxiality, in Sec. IV. 

The results of illustrative calculations are shown in 
Fig. 2, for X = 50, 100, and 200. The low-lying states 
form quasi-bands which may be roughly identified as a 
ground-state rotational band {K = 0), 7 vibrational ex- 
citation {K = 2), and two-phonon 7 excitations {K = 4 
and 0), denoted by 774 and 770. 

The stiffness of the potential around 7 = 0° simulta- 
neously determines both the 7-vibrational cnc^rgy scale 
]increasing from Fig. 2(a) to Fig. 2(c)] and also how well 
confined the wave function is with respect to 7, as seen in 
the corresponding approach to an ideal rotational spec- 
trum. Thus, within the framework of the Bohr Hamilto- 
nian, the 7 band energy more specifically, the energy 
ratio i?(2+)/i?(2^), or separation of vibrational and ro- 
tational energy scales — and the 7 softness of the wave 
function are inextricably linked. 

As a starting point, it may be observed that for x = 
the potential is strictly 7-independent, and the spectrum 
therefore follows an S0(5) multiplet structure [30, 31]. 
Successive multiplets consist of angular momenta 0, 2, 
4-2, 6-4-3-0, . . ., for w = 0, 1, 2, 3, . . ., respectively, with 
multiplet energies oc v{v + 3), as depicted in Fig. 3(a). 
The system is simply a Wilets-Jean [30] or SO (6) [32] ro- 
tor, but without (3 excitations (see also Ref. [14]). Then, 
as 7 confinement is introduced, the familiar rotational 
band structure begins to emerge. An intermediate spec- 
trum, obtained for x = 20, is shown in Fig. 3(b). 

For X = 50 [Fig. 2(a)], rotational quasi-bands are well- 
developed, and E{2j^)/E{2^) ss 10, as appropriate to, 
e.g., the well-deformed rare earth nuclei. However, it 
is seen from the potential plot in Fig. 2(a) that the 7 
confinement for this value of x is still weak. The range 
of energetically accessible 7 values increases significantly 
for successive phonon excitations, such that confinement 
is almost nonexistent at the energy of the two-phonon 
excitation. 

Dynamical 7 deformation consequently plays a major 
role in the calculated structure, through its interaction 
with the rotational dynamics. This is reflected in sig- 
nificant deviations from ideal rotational behavior in the 
spectroscopic predictions. 

Most noticeably, on inspection of Fig. 2(a), level ener- 
gies within the 7 quasi-band follow a gently 7-soft stag- 
gering pattern [2(34) (56) . . .]. This staggering is reminis- 
cent of the S0(5) level degeneracies obtained for x = 0, 
and it disappears as the 7 stiffness increases [Fig. 2(b,c)]. 
The deviations from rotational energy spacings are even 
more pronounced for the calculated two-phonon bands. 



4 




456789 10 i 456789 10 i 456789 10 ^ 



FIG. 2: Level schemes for the angular Hamiltonian (5) with ^ = 0, for (a) x ~ 50, (b) x ~ lOOj ''■nd (c) x ~ 200. Rotational 
1/(1/ + 1) energies for the yrast band are indicated by the dots. The potential V(7) is shown in the inset, with the ground, 
quasi-7, and quasi-77 bandhead energies indicated, (d-f) Staggering of level energies within the quasi-7 band, as measured by 
the energy second difference S{L). 
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FIG. 3: Level sehemes for the angular Harniltonian (5) with 
^ = 0, for (a) the 7-independent limit x — S'lid (b) x = 20, 
with levels arranged anticipating the quasiband structure of 
Fig. 2. Rotational L{L + 1) energies for the yrast band are 
indicated by the dots. The potential 1^(7) for x = 20 is shown 
in the inset, with the ground, quasi-7, and quasi-77 bandhead 
energies indicated. 

Note especially the near doubling of the rotational en- 
ergy spacing scale for the two-phonon bands, relative to 
the ground state band, for x = 50 [Fig. 2(a)]. 

The deviations from rotational energy spacings within 
the 7 band may be seen most clearly from plots of the 
level energy second difference S'(L) = \\E[L)—E{L—\)] — 
[E{L - 1) - E{L - 2)]\/E{2+), as shown in Fig. 2(d f). 
For an ideal rotational band with L{L + 1) energy spac- 
ings, the curve is flat, with S{L) = 1/3. Alternatively, 
7-soft staggering is manifest in minima at even L. As sur- 
veyed in Ref. [33], the observed level energies within the 
7 bands of most transitional and rotational nuclei yield 
S{L) plots which are either gently 7-soft or near constant 
(« 1/3). A few transitional nuclei {e.g., ^-^Sm, ^^''Gd, 
or ^^^Er) exhibit a degree of staggering comparable to 
that foimd for x = 50 (see also Rcfs. [17, 34]). However, 
most rare earth rotational nuclei (see Fig. 3 of Ref. [33]) 
more clearly follow a^a. L{L + 1) energy spacing within 
the 7 band. There is thus an apparent disagreement be- 
tween the degree of dynamical 7 softness expected in the 
Bohr picture given E{2^)/E{2f) w 10, and the observed 
structure in nuclei, at least if we assume the basic Hamil- 
tonian (5). 

Within the ground state band, the Hamiltonian (5) 
is found to yield relative energies [i.e., E{L^)/E{2^)] 
which fall below the L{L + 1) expectation for an adia- 
batic rotor. The ideal rotational energies are indicated, 
for comparison, by the dots in Fig. 2(a-c). The deviation 
from L{L + 1) spacing within the ground state band de- 
creases, as would be expected, for increasing 7 stiffness. 
The effect has already been noted in the context of a full 
^ and 7 calculation with the ACM in Ref. [13] (see Fig. 5 
of that reference). Such a deviation would traditionally 
be characterized as "centrifugal stretching" , based on an 



the interpretation in which the deformation increases, 
and thus the rotational moments increase, with increas- 
ing angular momentum. However, here the effect is seen 
to arise purely from the interaction of 7 and rotational 
degrees of freedom, for a system in which "stretching" in 
the /? degree of freedom is strictly impossible. 



B. Evolution of observables 

The evolution of the numerical predictions, with in- 
creasing 7 stiffness, is examined more quantitatively and 
systematically in Fig. 4. Both the energy spectrum 
[Fig. 4(lcft)] and electromagnetic (specifically, electric 
quadrupole) moments and transition matrix elements 
[Fig. 4(right)] are shown, as functions of x- 

The onset and evolution of rotational band structure, 
as 7 confinement is introduced, may be traced in the 
full energy spectrum [Fig. 4(a)]. Note especially the cor- 
relation between the 7 band energy [Fig. 4(a)] and the 
groimd state band energy ratio E{4^)/E{2f) [Fig. 4(b)], 
which varies from 2.5 for 7-independent rotation to 3.33 
for rigid axial rotation. This ratio is commonly taken as 
an indicator of rotational adiabaticity. For the present 
restricted problem, adiabaticity represents separation of 
the 7 and rotational degrees of freedom, but in general 
for the Bohr Hamiltonian the quantitative details will 
also be affected by the /3 degree of freedom. The evolu- 
tion of multiphonon band energies can also be followed 
in Fig. 4. These begin anharmonically low, at less than 
twice the 7 band energy — for x = 50, an estimate based 
on low-lying band members gives E^f-y^i/E-y w 1.7 and 
E^jfi/Ej ss 1.9 — but approach harmonicity as x in- 
creases. The relative energies of the bands may also be 
seen in Fig. 2(a-c). 

The evolution of electromagnetic properties is traced 
for representative quadrupole moments and transition 
strengths in Fig. 4(right). In the 7-independent limit, 
the wave functions are simply the SO (5) D SO (3) spher- 
ical harmonics themselves, and electromagnetic matrix 
elements are governed by S0(5) selection rules and re- 
lated by S0(5) D 80(3) Clebsch-Gordan coefficients. On 
the other hand, in the limit of large 7 stiffness, electro- 
magnetic matrix elements are expected to approach the 
Alaga rule ratios [2, 35] of the adiabatic axial rotor, given 
by ordinary angular momentum Clebsch-Gordan coeffi- 
cients. 

The electric quadrupole moments Q(2f) and Q(2+) 
are shown in Fig. 4(c). All quadrupole moments vanish in 
the 7-independent limit, by a selection rule arising from 
a parity quantum number defined in the five-dimensional 
space of the Bohr coordinates (M^-parity) [13, 21, 36]. In 
the rotational limit, these quadrupole moments are ex- 
pected to approach values of ±8-\/7r/7 ~ ±2.03, negative 
for the ground state band {K = 0) and positive for the 7 
band {K = 2), expressed relative to B{E2; 2+ 0^")^/^ 
These values are rapidly attained, by x ^ 25. 

For harmonic 7 vibration, the 7—^5, 774 — > 7, and 
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FIG. 4: Evolution of spectroscopic properties with 7 stifFness, for the angular Hamiltonian (5) with ^ = 0. Quantities shown 
are (a) excitation energies of low-lying levels, normalized to E{2f), (b) the energy ratio E{4:f)/E{2f), specifically, (c) electric 
quadrupole moments of the ground state band and 7 band 2+ members, (d) electric quadrupole reduced transition probabilities 
for one-phonon transitions between the ground, 7, and two-phonon 7 (/sT = and 4) bandhead states, and (e) reduced 
transition probabilities for the transitions depopulating the 2+ bandhead state. All electromagnetic quantities axe normalized 
to B{E2;2+ Of) = 1- 



770 ^ 7 iiitcirband intrinsic matrix elements {f\M'\i) [2] 
are expected to be in the proportion 1 : y/2 : 1 [18]. 
The overall normalization of these intrinsic matrix ele- 
ments, i.e., the 7 — > .9 strength, decreases with increasing 
7 stiffness [4]. For the transitions among the bandhead 
states, in particular, these intrinsic matrix element ra- 
tios correspond to B{E2;2+ 0+), B(i;2;4+^4 ^ 2+), 
and B{E2;0^ Q — > 2+) strengths in the proportion 
1 : 2.8 : 5. The approach to harmonic values is seen 
in Fig. 4(d). Simply from considering these transitions, 
harmonic behavior would appear to set in very gradually 
for X ^ 50. However, a more comprehensive consider- 
ation of the electromagnetic transition strengths, which 
leads to some modification of this conclusion, is provided 
by the Mikhailov analysis in Sec. HID. The branching 
ratios for electric quadrupole transitions between bands 
likewise approach the Alaga rule ratios. For the tran- 
sitions from the 2+ bandhead to the ground state band 



members [Fig. 4(e)], for instance, the adiabatic rotor has 
B{E2; 2+ ^ 0+), B{E2; 2+ ^ 2+) and B{E2; 2+ ^ 4+) 
strengths in the proportion 0.4 : 0.57 : 0.029. 



C. Effective 7 deformation 

Although we have so far examined 7 softness indirectly, 
through its spectroscopic signatures, the wave function 
■0(7, ri) is directly accessible for the eigenstates calcu- 
lated in the diagonalization of the Bohr Hamiltonian, and 
thus the deviation of 7 from 0° can be considered directly. 
The simplest measure is provided by an effective 7 value 
7, defined by 

cos 37 = (cos 37) . (10) 

The matrix elements of cos 37 in the S0(5) D S0(3) 
spherical harmonic basis are already available, as noted 
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FIG. 5: Evolution of the effective values 7 with respect to 
stiffness parameter Xi for the angular Hamiltonian (5) with 
5 = 0. Values are shown for ground state, 7, 774, and 770 
quasi-band members with L < 4. 



in Sec. II, so this expectation value may readily be 
calculated. The definition (10) is consistent with the 
quadrupole shape invariant approach [37, 38], in which 
an eflFective 7 for the full 7, fl) coordinate space is 
defined by cos37eff = (/J^ cos37)/(^2j,3/2 [39_4i]. 



The evolution of 7 for the ground state, 7, and 77 
band members (for L < 4) is shown in Fig. 5. In the 
X = (7-independent) limit, (cos 37) = by the M^- 
parity selection rule, and thus 7 = 30° for all states. 
As X increases past x ~ 50, it is seen that the 7 values 
for the members of each band cluster and decrease with 
increasing x- The 7 value jumps substantially between 
bands, increasing from ground to 7 to 77 bands, indeed, 
as expected for successive phonon excitations. 

The situation for "axial rotor" nuclei within the Bohr 
Hamiltonian framework is very much contrary to the clas- 
sic but schematic characterization of such nuclei as hav- 
ing "7 « 0°" , which may be more concretely interpreted 
as 7 <C 30°. Recall that the 7-band excitation energies 
matching the experimental values for rotor nuclei arc ob- 
tained for X ~ 50. For this stiffness, the ground state 
band members have 7 « 15°, and the 7 band members 
have 7 w 23°. These large 7 values are consistent with 
the large range of energetically accessible 7 values for 
these states [Fig. 2(a,inset)]. The full probability distri- 
bution with respect to the 7 coordinate is considered in 
Sec. V. 



D. Intrinsic matrix elements 

A more comprehensive and meaningful examination of 

electromagnetic transition strengths is realized by consid- 
ering the inter band transitions in aggregate, according to 
the Mikhailov mixing formalism [42] . Within this frame- 
work, all transition amplitudes are expressed in terms of a 
single intrinsic electromagnetic matrix element and single 
mixing parameter between each pair of bands. The am- 
plitudes are expected to fall on a straight line on an ap- 
propriate (Mikhailov) plot of {K2J2\\-M\\KiJi) or, com- 
monly, B{E2f/'^ vs. J2(J2 + 1) - ^"1(^1 + !)■ The intrin- 
sic matrix elements and mixing parameter are identified 
from the slope and intercept. 

Specifically, for interband transitions with AK = 2, 
the leading-order band mixing relation for E2 reduced 
matrix elements is [2, (4-210)] 

{K2J2\\M\\KiJi) = ai(2Ji + l)i/2(Jiifi22|J2X2) 

X [Ml + M2[J2(J2 + 1) - MJi + 1)]] , (11) 

where it is assumed that K2 — Ki+2, and where cti = y/2 
if Ki = or CTi = 1 otherwise. The parameters in this 
expression are related to the intrinsic matrix element 
{K2\M'\Ki), mixing matrix element {K2\s+2\Ki), and 
intrinsic quadrupole moment Qq by Mi = {K2\A4'\Ki) — 
A{Ki + 1)M2 and M2 = [15/{8TT)]^/^eQo{K2\e+2\Ki) [2, 
(4-211)]. The intrinsic matrix element may thus be ex- 
tracted from the slope and intercept as 



{K2\M'\Ki) = Ml + 4{Ki + 1)M2. 



(12) 



More specific expressions for K- decreasing and K- 
increasing transitions, in terms of B[E2) reduced tran- 
sition probabilities, are given in Appendix B. 

The interband quadrupole transition strengths for the 
Bohr Hamiltonian calculations of Sec. Ill A are shown in 
Fig. 6 in Mikhailov form. They are plotted as ^(£2)^2 
vs. Lf{Lf + 1) — Li{Li + I), for transitions between 
states with L < 6. For the most part, the transition 
amplitudes do indeed follow an essentially linear pat- 
tern, and it is therefore meaningful to extract effective 
intrinsic matrix elements, well as mixing parameters, 
from the Mikhailov analysis. (The Mikhailov formalism 
has been applied to extract effective intrinsic matrix el- 
ements from the interacting boson model [43] , in a simi- 
lar fashion, in Rcfs. [44, 45].) However, deviations from 
a linear relation are significant for transitions involving 
the two-phonon quasi-bands for x = 50 [Fig. 6(left)], 
as might be expected from the substantial 7-softness 
and deviations from rotational energy spacings already 
noted for these bands. The resulting intrinsic matrix 
elements for the 7 — >■ 774 — >■ 7, and 770 — > 7 tran- 
sitions, obtained from (B2) and (B4), are listed in Ta- 
ble I, together with the dimensionless mixing parame- 
ter a = IM2/M1I (see Appendix B). The normalization 
of the electric quadrupole operator M.{E2) is arbitrary 
in the present analysis. To provide a scale for compar- 
ison with experiment, the intrinsic matrix elements in 
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FIG. 6: Interband transition amplitudes B{E2)^^^ , from the 7 quasi-band to the ground state band (top), from the 774 quasi- 
band to the 7 quasi-band (middle), and from the 770 quasi-band to the 7 quasi-band (bottom), for Mikhailov analysis. Plots 
are included for the calculations of Fig. 2, with x = (l^'ft). X = 100 (middle), and x ~ 200 (right) and ^ = 0. The values 
shown are for transitions between levels with i < 6, normalized to B{E2; 2^ — > Of) = 1. 



TABLE I: Electric quadrupole interband intrinsic matrix elements {f\M'\i) and mixing parameters a, for different 7 stiffnesses, 
as extracted from the Mikhailov analyses of Fig. 6. Ratios, as indicators of anharmonicity, are tabulated in the final two 
columns. The values for an adiabatic rotor with harmonic 7 vibration [18] are included for comparison. The values for the 
intrinsic matrix elements are normalized to B{E2- 2f — >■ Oj") = 1. 
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" The 77 —^ 7 intrinsic matrix elements for x = 50 can only be crudely approximated, since the Mikhailov plot yields values 
which are not strongly linear [Fig. 6(d,g)]. The estimated parameters used in the analysis are Mi fs 0.9 for 774 7 and 

Ml « 0.4 for 770 — >■ 7. 
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Tabic I are given relative to the square root of the in- 
band B{E2;2+ 0+). 

For harmonic 7 vibration, the ratios of the 77 — >■ 7 in- 
trinsic matrix elements to the 7 — > g intrinsic matrix el- 
ement are expected to be (7|A^'|774)/(5|A^'|7) = -\/2 « 
1.41 and (7|A^'|77o)/(5|A^'l7) = 1, according to the pro- 
portion noted in Sec. IIIB. For comparison, ratios of the 
intrinsic matrix elements extracted from the Bohr Hamil- 
tonian numerical calculations arc given in the last two 
columns of Table I. Note the rapid quantitative approach 
of these calculated ratios to the expected harmonic val- 
ues. Even the 77 — > 7 transitions for the soft x = 50 
case are essentially consistent with harmonic ratios, to 
the extent that slope and intercept parameters can mean- 
ingfully be extracted in this instance [Fig. 6(d,g)]. For 
X = 200, harmonic values are obtained to within ^ 1%. 

The bandmixing, indicated by the Mikhailov plot 
slopes, is substantial in all the cases considered in Ta- 
ble I. The harmonicity of the intrinsic matrix elements 
is therefore not apparent simply from the plot inter- 
cepts buy only after the leading-order bandmixing cor- 
rections (B2) and (B4) are taken into account. For ex- 
ample, even for the most adiabatic case, x ~ 200, the 
7 — >■ g [Fig. 6(c)]and 774 — > 7 [Fig. 6(f)] Mikhailov plots 
both have slope parameters a w 0.012, resulting in a 5% 
adjustment to the 7 — >■ g intrinsic matrix element and a 
14% adjustment to the 774 — > 7 intrinsic matrix element. 

In summary, although the strengths of the individ- 
ual interband transitions only approach the limit of an 
adiabatic rotor (and, more spccificially, harmonic vibra- 
tion) gradually, as observed from Fig. 4(d), this devia- 
tion is quantitatively well-described in terms of a rapid 
approach to harmonic values of the interband intrin- 
sic matrix elements, but with the individual transition 
strengths modified by leading-order AK = 2 bandmix- 
ing (11). The strength of this mixing then gradually 
decreases with increasing 7 stiffness. 



IV. ONSET OF RIGID TRIAXIALITY 

The excitation spectrum may be expected to change 
dramatically with the onset of rigid triaxiality. The Bohr 
Hamiltonian predictions ultimately approach a 7 = 30° 
Davydov rotor spectrum [8] for confinement by a suffi- 
ciently stiff cos^ 87 potential [13]. However, the initial 
onset of triaxiality is reflected in much more subtle de- 
viations from the characteristics of an axially symmetric 
rotor. The difference between axial and triaxial minima 
in the potential is obscured by the substantial dynami- 
cal fluctuations in 7 present in both cases. As noted in 
Sec. II, the onset of triaxiality may be investigated by 
considering the introduction of a cos^37 contribution, 
i.e., nonzero ^, in the Hamiltonian (5). 

The results of calculations for two representative po- 
tentials are shown in Fig. 7: the soft or "critical" ax- 
ial minimum (^ = 0.5) [Fig. 7(b)] and a weakly triax- 
ial minimum (^ = 0.8) [Fig. 7(c)]. For each of these 



calculations, the potential depth, or x, is chosen to 
give E{2^)/E{2^) « 10, again appropriate to the well- 
deformed rare earth nuclei. The comparable axial rotor 
calculation with the same 7 band energy, i.e., x = 50, is 
shown again as a baseline for comparison [Fig. 7(a)]. 

In Fig. 7, the 7-phonon quasiband structure is seen to 
remain intact. Our concern is therefore with the prin- 
cipal spectroscopic properties of these bands — excita- 
tion energies of the bands, deviations from rotational en- 
ergy spacing within the bands, and electric quadrupole 
intrinsic matrix elements. The two-phonon energy anhar- 
monicities evolve from slightly negative {E^^/E^ < 2) for 
^ = [Fig. 7(a)] to positive {E^-y/E^ > 2) [Fig. 7(b,c)] 
with the introduction of triaxial tendencies. The anhar- 
monicity of the 770 band rises more rapidly than that of 
the 774 quasi-band. Qualitatively, this is consistent with 
evolution towards a 7-stiff, adiabatic triaxial rotor [18], 
for which the K = 4 quasi-band is a triaxial rotational 
excitation and the K = quasi-band is a 7 vibrational 
excitation. 

The level energies within the 7 band progress, with 
increasing from 7-soft staggering [2(34)(56) . . .] to 
the reverse pattern associated with triaxial rotation 
[(23)(45) . . .] [8]. As in Sec. Ill A, the staggering may 
be seen most immediately from plots of the second dif- 
ference S{L) [Fig. 7(d-f)], which has minima at even L 
for 7-soft staggering or at odd L for triaxial staggering. 

The "centrifugal stretching" phenomenon in the yrast 
band, i.e., reduction of £^(L^)/i?(2^) relative to L(L-fl) 
spacing, persists [Fig. 7(b,c)] at about the same level as 
for X = 50. However, the growth in rotational constant 
(and general deviation from rotational behavior) for the 
excited, especially 77, bands is tamed relative to the ax- 
ial calculation. This may be at least qualitatively un- 
derstood by comparing the potential plots in Fig. 7(a- 
c, insets). The axial calculation of Fig. 7(a), as noted in 
Sec. HI A, provides only weak confinement at the 77 band 
energies (7 < 40°). Although the nominally "softer" cal- 
culation of Fig. 7(b) does provide weaker confinement, 
compared to this axial calculation, at the ground state en- 
ergy, it actually provides stiffer confinement, to a smaller 
range of 7 values (7 < 30°), at the 77 band energies. 
[This effect may be more properly considered a reflec- 
tion of the steep rise in the cos^ 87 term used to create 
the triaxial confinement than an intrinsic property of the 
onset of triaxiality per sc. There is no inherent calcu- 
lational reason not to consider a potential with, for in- 
stance, a triaxial minimum located at the same position 
as in Fig. 7(c, inset) but a lower barrier at 7 = 60°.^] A 
similar observation may be made for the calculation of 
Fig. 7(c), which provides confinement to triaxial 7 at the 



Any potential V(-f) satisfying the basic requirements from the 
Bohr coordinate symmetries may be expanded in terms of the 
form cos" 87 (which is equivalent to Fourier decomposition in 
terms of the form cos Sn^y) and therefore may readily be accom- 
modated for calculations within the ACM. 
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FIG. 7: Level schemes for the angular Hamiltonian (5), for (a) ^ = with x = 50, (d) ^ = 0.5 with x = 100, and (e) f = 0.8 
with X = 500. Rotational L{L + 1) energies for the yrast band are indicated by the dots. The potential 1^(7) is shown in the 
inset, with the ground, quasi-7, and quasi-77 bandhead energies indicated, (d-f) Staggering of level energies within the quasi-7 
band, as measured by the energy second difference S{L). Figure adapted from Ref. [21]. 



TABLE II: Electric quadrupole interband intrinsic matrix elements {f\A4'\i) and mixing parameters, for different 7 potential 
shapes chosen to reproduce the onset of weak rigid triaxiality, as extracted from the Mikhailov analyses of Fig. 8. Ratios, as 
indicators of anharmonicity, arc tabulated in the final two columns. The Y(5) triaxial estimate [10] is included for comparison. 
The values for the intrinsic matrix elements are normalized to B{E2;2f — >■ Oj) = 1. The results in this table also serve 
to correct intrinsic matrix element values given previously in Table 1 of Ref. [26]. The roles of initial and final bands were 
interchanged, in that analysis, when extracting the slope parameter from (Bl) and (B3), resulting in the use of an incorrect 
sign for the bandmixing correction term in (B2) and (B4). 
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ground state energy, but simply provides (axial) confine- 
ment to 7 < 30° at the 77 band energies. 

For the weakly triaxial calculations considered here, 
the interband transition strengths continue to follow an 

essentially linear pattern on a Mikhailov plot, as expected 
for rotational bandmixing, as seen in Fig. 8. The 77 —> 7 



transitions, in fact, demonstrate better linear behavior 
[Fig. 8(e-f,h-i)] than for x = 50 [Fig. 8(d,g)]. Inter- 
band intrinsic matrix elements may therefore again be 
extracted from the Mikhailov analysis, as given in Ta- 
ble II. The intrinsic matrix element remains es- 
sentially constant, and equal to that for the axial x = 50 
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FIG. 8: Interband transition amplitudes B{E2)^^'^ , from the 7 quasi-band to the ground state band (top), from the 774 quasi- 
band to the 7 quasi-band (middle), and from the 770 quasi-band to the 7 quasi-band (bottom), for Mikhailov analysis. Plots 
are included for the calculations of Fig. 7, with ^ = (x = 50) (left), ^ = 0.5 (x = 100) (middle), and ^ = 0.8 (x = 500) (right). 
The values shown are for transitions between levels with L < 6, normalized to B{E2;2f Of) = 1. Figure panels (a-f) 
adapted from Ref. [21]. 



calculation, but the 774 — >■ 7, and 770 — > 7 intrinsic 
matrix elements decrease substantially compared to the 
harmonic 7- vibrational values. 



Such a reduction of the 77 — )■ 7 intrinsic matrix ele- 
ments, relative to the harmonic values, has already been 
proposed [10] on relatively simple grounds. Supposing 
an adiabatic separation of rotation from vibration, and 
furthermore imposing a small-7 approximation, yields a 
one-dimensional Schrodinger equation problem in 7. In 
Ref. [10], a square well is then adopted for 1^(7) to simu- 
late the onset of triaxiality. This yields the Y(5) estimate 
shown for comparison in Table 11. 



V. WAVE FUNCTION PROBABILITY 
DISTRIBUTIONS 

In the limit of adiabatic separation of the 7 and rota- 
tional degrees of freedom, the wave functions of all mem- 
bers of a band would be given by 



(13) 



where the function Fk (7) would be identical for all states 
within the same band, independent of L. The band is 
characterized by intrinsic angular momentum projection 
K. This may be contrasted to the general situation (6), 
in which all even K with < K < L {or 2 < K < L ioT L 
odd) can contribute, and the coefficients Fk (7) need not 
be directly related for different states. The breaking of 
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FIG. 9: Probability distributions with respect to 7 for low-lying quasi-band members in caleulations with Hamiltonian (5), for 
the ax-iid cases x = 50 (left) and x = 200 (middle), both with ^ = 0, and for the weakly tria:xial case { = 0.8 with x = 500 
(right). Probability distributions are shown for members of the ground state, 7, 774, and 770 quasi-bands (top to bottom, 
respectively), with L < 6. 



adiabaticity has already been seen to have spectroscopic 
consequences (Sees. Ill and IV). Here we shall more di- 
rectly inspect the wave functions themselves, through the 
probability distributions. Specifically, we examine the 
probability distribution P(7), with respect the 7 coordi- 
nate, after integration over Euler angles, and the proba- 
bility decomposition Pk , with respect to the K quantum 
number for the Euler angle (rotational) dependence, af- 
ter integration over 7. The calculational details are given 
in Appendix C. 

First, considering -P(7), results are given in Fig. 9 
for the softest axial calculation of Sec. Ill (x = 50) 
[Fig. 9(lcft)], the stiffest axial calculation of Sec. Ill 
(x — 200) [Fig. 9(middlc)], and the weakly triaxial cal- 
culation of Sec. IV (^ = 0.8 with x = 500) [Fig. 9(right)]. 
Successive panels (top to bottom) show the ^(7) distri- 
butions for the ground, 7, 774, and 770 band members, 
respectively, with L < 6. All the ^(7) vanish at 7 = 0° 
and 7 = 60°, due to the volume element for the Bohr 
coordinates (see Appendix C). 

The basic features seen in Fig. 9 may be qualitatively 



understood in terms of the small-7 limit of (5), which re- 
duces {e.g., Ref. [10]) to a two-dimensional harmonic os- 
cillator problem, with two-dimensional angular momen- 
tum TO = K/2 and with 7 as the "radial" variable. The 
K = 2n-y (or to = n^) bands, i.e., the ground, 7, and 774 
bands, have probability distributions which are nodeless. 
These move towards higher 7 with increasing phonon 
number n-y [Fig. 9(a,d,g) or Fig. 9(b,e,h)]. The cen- 
ters of the probability distributions are at substantially 
nonzero 7 values, in the 10°-30° range, but move towards 
smaller 7 for larger stiffess [compare Fig. 9 (left) with 
Fig. 9(middle)]. All these properties are as anticipated 
from the 7 values in Fig. 5. For the 770 band, which is 
characterized by if = 2{n^ — 2) (or m = — 2), the 
probability distributions have a single node [Fig. 9(j,k)]. 

Adiabatic separation (13) implies identical P(7) dis- 
tributions for all members of the same band. Indeed, 
the -P(7) curves are virtually indistinguishable between 
band members for the examples in Fig. 9. The exceptions 
are, once again, the 77 bands in the x = 50 calculation 
[Fig. 9(g,j)]. There is some slight displacement between 
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FIG. 10: The K content of low- lying quasi-band members in calculations with Hamiltonian (5), for the axial cases x = 50 (left) 
and X = 200 (middle), both with ^ = 0, and for the weakly triaxial case ^ = 0.8 with x = 500 (right). Probabilities Pk for 
= (solid curve), K = 2 (dashed curve), and K = ^ (dotted curve) are shown for members of the ground state, 7, 774, and 
770 quasi-bands (top to bottom, respectively), with L < 10. 



the curves for the different members of the ground or 7 
bands in this calculation as well. The breaking of adia- 
baticity is also apparent for the 770 band members with 
L > 0, from the disappearance of the node in P{j), which 
indicates that multiple K values must contribute to the 
wave function.^ 

It is interesting to note the qualitative differences of 



^ When only one K term contributes to (6), as in (13), a zero- 
crossing in Fk(j) necessarily yields a zero- valued minimum in 
Pi'j)- If, instead, the minimum is washed out, it may be con- 
cluded that multiple K terms are contributing in (C2), such that 
these terms do not simultaneously have nodes at the same 7 
value. 



the more triaxial calculation [Fig. 9(right)] from the ax- 
ial calculations [Fig. 9(left, middle)]. The P{j) for the 
ground, 7, and 774 bands {i.e., those with nodeless dis- 
tributions) [Fig. 9(c,f,i)] are peaked at 7 values roughly 
comparable to those for the x = 50 "axial" calculation 
[Fig. 9(a,d,g)] (recall that the parameters were chosen so 
that these calculations share the same 7 band energy) but 
are more sharply peaked. The 770 distribution [Fig. 9(1)] 
shows a marked enhacement of the peak at small (axial) 
7. This may seem counterintuitive for a "triaxial" calcu- 
lation, but, as already remarked in Sec. IV, the triaxial 
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confinement is limited to the ground state band energy.^ 
In interpreting the -P(7) distributions as indicators of 
adiabaticity, it should be noted that, although adiabatic 
separation implies identical Pi"/) distributions, the con- 
verse is not strictly true. Adiabaticity might be violated, 
and several K values might contribute in (6), but the 
various Fk (7) for the different band members may be re- 
lated such that, nonetheless, the same ^(7) distributions 
are obtained after integration over Euler angles. There- 
fore, these distributions can only be conclusively taken 
to indicate adiabaticity if it is also known that only one 
K value contributes significantly. 

The contributions of different K values in each of the 
bands (ground, 7, 774, and 770) are shown in Fig. 10, 
for each band member with L < 10. For the calculations 
in Fig. 10, the bandhead states have essentially pure K. 
The largest admixture in a bandhead state is ~ 3% for 
the 774 bandhead in the x = 50 calculation, but the 
bandhead K admixtures in the other calculations are all 
< 10~^. (Note that the 770 bandhead, as an L = 
state, trivially has pure K = 0.) The admixtures in- 
crease with L within each band. Again, the extremes are 
in the 77 bands for x = 50, where the admixtures ac- 
count for approximately half the probability at L = 10 
[Fig. 10(g,j)]. In contrast, for the weakly triaxial calcu- 
lation [Fig. lO(right)], the K admixtures in the 77 bands 
are actually slightly smaller than for the ground state 
band. Indeed, they closely match the K admixtures of 
the corresponding bands in the stiff axial x = 200 calcula- 
tion [Fig. lO(middle)]. This observation is consistent with 
the characterization of these bands as relatively "good" 
axial rotational bands, as suggested spectroscopically in 
Sec. IV. 



VI. CONCLUSION 

The possibility of exact diagonalization of the Bohr 
Hamiltonian for essentially arbitrary /3 and 7 stiffness 
opens the door for direct comparison of the Bohr Hamil- 
tonian predictions with experiment throughout the range 
of possible dynamics for the nuclear quadrupole degree 
of freedom. At a phenomenological level, this permits 
meaningful tests of the Bohr Hamiltonian for general 
rotor- vibrator nuclei. 

For instance, in the past, interpretation of rotational 
"phonon" states, although nominally within the Bohr 
description, has largely been at a schematic level {e.g., 
Refs. [44, 46-49]): adiabatic separation of the rotational 



Moreover, under adiabatic separation, the 7 wave function for 
the excited K = hand must be orthogonal to the ground state 
band wave function. Since this distribution has moved to larger 
7 values, the redistribution in probability to the smaller-7 peak 
for the excited band can be understood from orthogonality con- 
straints, following arguments similar to those applied in Ref. [29] 
for prolate-oblate coexistence. 



and vibrational degrees of freedom is assumed, the j3 and 
7 excitations are taken to be harmonic, and phonon se- 
lection rules are assumed for electric quadrupole transi- 
tions. These predictions are then adjusted by the leading- 
order spin-dependent bandmixing relation, but with ad 
hoc mixing parameters. Here, instead, we explore exact 
predictions of the Bohr Hamiltonian, both for axial and 
weakly triaxial confinement. 

The present analysis, which has been restricted to the 
7 and rotational degrees of freedom, provides a starting 
point for understanding the full dynamics involving all 
five Bohr coordinate degrees of freedom, i.e., consider- 
ing coupling with the /3 degree of freedom as well. Many 
of the qualitative properties of the present solution may 
be expected to carry over (see, e.g., Fig. 4 of Ref. [26]). 
However, the introduction of /? softness may generally be 
expected to quantitatively alter the results, for instance, 
further attenuating the rotational character of the bands 
[e.g., reducing the ratio £'(4j^)/_E(2j^)]. Moreover, in the 
case of near degeneracy of the 7 phonon or multiphonon 
bands with bands involving (3 excitations, bandmixing 
can substantially alter the results. Therefore, detailed 
comparison with experiment should be made in the con- 
text of a full treatment incorporating /3 softness. 

Microscopic descriptions of nuclear collectivity rely 
upon a reduction of the many-body problem to one in- 
volving effective collective degrees of freedom. Mean- 
field approaches to deriving the quadrupole collective 
dynamics (reviewed in, e.g., Refs. [27, 50, 51]) yield a 
Hamiltonian involving a much more general, coordinate- 
dependent form for the kinetic energy operator than the 
conventional but schematic Laplacian form considered 
in (1). The resulting generalized Bohr Hamiltonian [27] 
may be represented in terms of coordinate-dependent mo- 
ments of inertia. It should be noted that the ACM can 
readily accommodate Hamiltonians involving much more 
general differential operators [16] in the (3 and angular 
variables than the simple Laplacian form. For instance, 
scalar-coupled products of the quadrupole momentum 
tensor p and coordinate tensor q constitute an impor- 
tant special case considered in the geometric collective 
model [4, 6]. The Bohr kinetic energy is obtained as the 
lowest-order term {pxpY'^\ and attention in phenomeno- 
logical studies has largely been limited to the next term 
{p X q X pY^\ These and higher-order terms in the co- 
ordinate dependence may be combined to recover much 
or all of the flexibility of the generalized Bohr Hamilto- 
nian [52]. 

Even further generalizations may be required. For in- 
stance, the Sp(3,K) symplectic shell model framework 
gives rise to a collective model in which the generalized 
Bohr Hamiltonian must be augmented with a vorticity 
degree of freedom [53] . Since the collective model serves 
as the intermediate link between microscopic theories and 
spectroscopic predictions, it is essential to determine the 
limitations of the Bohr Hamiltonian and the nature of 
the modifications required such that its predictions can 
accurately describe the observed phenomena. 
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Appendix A: Restriction to angular coordinates 

In this appendix, the reduction of the full Bohr Hamil- 
tonian (1) to an angular Hamiltonian (4) is briefly sum- 
marized. First, for convenience, let us simplify the Bohr 
Hamiltonian to its equivalent dimensionless form 



H = - 



A! 

13' 



where 



(Al) 



(A2) 



This appendix adapts the leading-order AK = 2 band- 
mixing relations (11) and (12) to the form required for 
the analysis of Figs. 6 and 8. For K- decreasing transi- 
tions {e.g., 'Y ^ g and 774 7), in terms of B{E2) 
reduced transition probabilities, 

B{E2-KiJi ^ KfJf) = af{JiKi2 - 2\JfKff 

X Mf[l + a[Jf{Jf + l)-Ji{Ji + l)]]\ (Bl) 

with normalized positive slope parameter a ~ — M2/M1. 
Thus, the intrinsic matrix element is extracted as 

{Kf\M'\Ki)=Mi[l-4{Kf + l)a]. (B2) 

Similarly, for K-increasing transitions {e.g., 770 7), 

B{E2; KiJi ^ KfJf) = aUJiKi2 + 2\JfKff 

X M^[l + a[Jf{Jf + 1) - Ji{Ji + l)]]^ (B3) 

where now the positive slope parameter is a = +M2/M1, 
and thus the intrinsic matrix clement is extracted as 



{Kf\M'\Ki) = Mi[l + A{Ki + l)a]. 



(B4) 



by rescaling H {2B/h'')H and V {2B/h')V. The 
two routes to obtaining an angular Hamiltonian indicated 
in Sec. HA proceed more precisely as follows: 

(1) Schematically, rigid (3 deformation (/3 « /?□) is ob- 
tained if the nuclear wave function ^'(/3,7,f2) is highly 
localized by a stiff potential with respect to p. For speci- 
ficity, consider V{P,j) = u{(3) + v{"i). Then H Ki Hfj + 
Po'H^vt, where Hp = -K+u{(3) and if^o = k'+Plv{-i). 
The separated eigenfunctions \l/(/3,7, O) = /(/3)V'(7, fi) 
satisfy Hpf{f3) = efjf{f3) and H^nil;{j,n) = ejQip{j,fl). 
Note that the angular problem is thus of the form (4), 
with 1^(7) = (3qv{'j). The total energy eigenvalues E, 
defined by iJ*(^,7,il) = £;*(/3,7, f2), are obtained ad- 
ditively as E = ep + P^'^s^n. Therefore, for fixed (3 
excitation {e.g., the ground state for the /3 problem), the 
eigenvalues of the angular problem directly give the en- 
ergy spectrum. These arguments apply only in the limit 
of stiff /3 confinement, and finite /3 softness may be ex- 
pected to lead to /3-7 coupling [17]. 

(2) Alternatively, for V{P,'y) = u{p) + v{-f)/(3', 
the Bohr Hamiltonian eigenproblem is exactly sepa- 
rarable [28]. In this case, H = Hp + where 
Hjj = -A + u(/3) and now H^q = + v{j). The 
separated eigenfunctions *(/?, 7, il) = f{l3)tp{j,V,) sat- 
isfy {Hp + f3-'e^n)f{P) = Ef{/3) and H^n^P{-f,n) = 
Sjn'>P{l,^)- Note that the angular problem is of the 
form (4) with ^(7) = v{j). The eigenvalue s^q from 
the angular problem now appears in the /3 equation as a 
"centrifugal" coefficient, i.e., multiplying /3~^. It there- 
fore enters indirectly into the total eigenvalue E, through 
the /3 eigenproblem, rather than directly giving the en- 
ergy spectrum. 



Appendix C: Wave function probability relations 

In this appendix, expressions are given for the prob- 
ability distribution P{j), with respect to the 7 coordi- 
nate, and the decomposition Pk, with respect to the K 
quantum number, for a wave function ^"(7) Note that 
the volume element for the coordinates (7, fl) is given by 
[sin 37I djdil,. 

The probability distribution ^(7) is obtained by inte- 
gration over Euler angles, as 



P(7) = |sin37| \i;{^,n)\'dn. 



(CI) 



and thus, in terms of the (real) coefficient functions 
Fx (7) appearing in (6), 



IStt 



(02) 



K=0 



The angular integration has been carried out using the 
orthogonality integral for the functions [54], which 

gives /4U'(f^)4M(f^)rf^ = [(167r2)/(2L + l)]5^a 
5k' k5m' M ; unless K = Q with L odd, in which case the 
integral vanishes [21]. For the eigenfunctions tpLiM{l, ^) 
obtained with respect to the S0(5) D S0(3) basis, the 
known quantities are the diagonalization coefficients OLij 
appearing in (9) and the functions FLiK{7) in the repre- 
sentation of the S0(5) D S0(3) spherical harmonics 



Um{1,^) = Yl PLiK{l)£!Klm^ 
even 



(03) 
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where we again use a counting index to label the spherical 
harmonics. In terms of these, ^ 

PlAi) = 2^^-^ Ism 37 1 

L 

X ^^aLijaLikFLjK{l)FLkK{l)- (C4) 

K=0 jk 
even 

The contribution of each K value to ^^(7, f^), integrated 
over 7, is 



167r2 



2L 



T / [^'^(7)]^ sin 37^7. (C5) 
^ Jo 



* The expression (C4) for P(7) is equivalent to (A7) of Ref. [17]. 
However, the normalization factors appearing in these expres- 



For the functions ipLiMi"/,^), represented by a^ij coef- 
ficients with respect to the S0(5) D S0(3) basis, these 
probabilities may be computed as 



O-Lik 



jk 



X / FLjKil)FLkKh)sm3jdj. (C6) 
Jo 



sions differ, due to the different normalization conventions de- 
fined for the Fx(7) coefficients in Ref. [17] and in the present 
work, which follows Ref. [21]. 
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